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Abstract 

An non-autonomous system is proposed to model the seasonal pat- 
tern of dengue fever. 

We found that an approximate threshold condition for infection 
persistence describes all possible behavior of the system. 

As far as we know, the kind of analysis here proposed is entirely 
new. No precise mathematical theorems are demonstrated but we give 
enough numerical evidence to support the conclusions. 
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1 Introduction 

In a previous paper pQ, in which we tried to understand the phenomenon 
of dengue overwintering, we discovered an interesting threshold condition 
that allows the complete qualitative understanding of the behavior of a non- 
autonomous system. 

To motivate the reader we briefly describe the phenomenon we studied in 

m 

In subtropical regions dengue fever, a mosquito transmitted disease, shows 
a resurgent pattern with yearly epidemics, which starts typically in the 
months characterized by heavy rains and heat, peaking some three or four 
months after the beginning of the rainy season. In the dry months the num- 
ber of cases drop essentially to zero due to the virtual disappearance of the 
vector. Since the infection reappears for some years in the same regions, it 
is natural to ask how the virus survives the dry season. 

In order to model those seasonal patterns of the disease, we proposed a 
non-autonomous system, described below. 

The model describes the dynamic of dengue in its three components of 
transmission, namely, human hosts, mosquitoes and their eggs (the latter 
includes the intermediate stages, like larvae and pupae). These populations, 
in turn, are divided into susceptible humans, denoted Sh, infected humans, 
Ih, recovered (and immune) humans, Rh, total humans, N H = Sh+Ir+Rh, 
susceptible mosquitoes, Sm, infected and latent mosquitoes, L M , infected and 
infectious mosquitoes, Im, non-infected eggs, Se, and infected eggs, Ie- 

The model's dynamics is described by the set of equations. 
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= -abI M $fc-VL H S H + r H N H (l-!j£) 

= abI M f^-( l 2 H + a H + lH )I H 

^ = lH I H -pL H R H 

^ff = Ps (cs ~ d s sin (27r/* + 7r)) S E 6 (c s - d s sin (27r/* + 7r)) 

dJ T = aSM^-e-^aSM(t-TT)^^-iM M L M 

pj (c/ — sin (27r/t + 7r)) i^fl (cj — sin (27r/t + 7r)) 

^ = [r M ^M + (1 - g) t m Im\ (l - ^i) - 

fJ'E^E ~ Ps {cs - ds sin {2itft + 7r)) S^fl (c e - d E sin (27r/£ + 7r)) 



^ = ^m/m (l - <Sgfel) - 

(c/ — <i/ sin (27r/t + 7r)) (c/ — dj sin (27r/t + 7r)) 

(1) 

Let us briefly describe some features of the model. 

We begin by describing the first three equations of the model. 

Susceptible humans grow at the rate thNh(1 — ^) — [i h Sh, where r# 
is the birth rate, fi H is the natural mortality and is the human carrying 
capacity. Note that we are assuming that close to the carrying capacity 
the human population growth is checked by a reduction in the birth rate. 
Alternatively the control of the population could be done by assuming an 
increase in the mortality rate, but the net result would be qualitatively the 
same. Those susceptible humans who acquire the infection do so at the rate 
ablMjf^, where a is the average daily biting rates of mosquitoes and b is 
the fraction of infective bites inflicted by infectious mosquitoes Im- Infected 
humans, In may either recover, with rate 7, or die from the disease, with 
rate (fi H + an)- Recovered humans remain so for the rest of their lives. 

The fourth, fifth and sixth equations represent the susceptible, latent and 
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infectious mosquitoes sub-populations, respectively. Susceptible mosquitoes 
varies in size with a time-dependent rate 

p s (c s - d s sin (2nft + tt)) S E 9 (c s - d s sin (2nft + tt)) - /i M S M - 

The term [i M is the natural mortality rate of mosquitoes. The term psSe is 
the fraction of eggs present at time t, and which survived the development 
through the intermediate stages (larvas and pupas). The time-dependent 
rate (q — d; L sin (2nft + vr)) 9 (q — d; L sin (2n ft + n)) (i — S, I) simulates the 
seasonal variation in mosquitoes production from eggs, assumed different for 
infected and susceptible eggs, for generality. By varying q and di ,(i — S, I), 
we can simulate the duration and severity of the winters (/ = 1/365 days -1 
and so it fixes one cycle per year). The Heaviside ^-function (a step function 
that is equal to zero when the argument is less than zero and one when the 
argument is greater or equal to zero) 6 (q — djsin (27r ft + 7r)) prevents the 
term 

(ci - di sin (2nft + n)) 9 (q - di sin (2nft + tt)) (i = S, I) 
from becoming negative. If q is smaller than di, then the winter is long 
and severe. On the other hand, if q is greater than di, then the winter is short 
and mild. Susceptible mosquitoes who acquire the infection do so at the rate 
uSmj^, where a is the average daily biting rates of mosquitoes, and became 
latent. A fraction of the latent mosquitoes survives the extrinsic incubation 
period with probability e~^ MTI and become infectious. Therefore, the rate 
of mosquitos becoming infectious per unit time is e~^ MTI aSuit — Tj) • 
The term 



represent vertical transmission, that is, the rate by which infected eggs be- 
come infectious adults. Infected mosquitoes die at the same rate fi M as the 
susceptible ones. 

The seventh and eighth equations represent the dynamics of susceptible 
and infected eggs, respectively. 

In the seventh equation, the term 



represent the birth rate of susceptible eggs born from susceptible mosquitoes 
with rate 



pi (cj — dj sin (2nft + tt)) Ie9 (cj — dj sin (2irft + tt)) 
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and from a fraction (1 — g) of infected mosquitoes, with rate 



The term tm 



1 — - g , E> ) represents the density-dependent rate of eggs 




birth rate. Once again we choose a density dependence on birth rather than 
on death. Alternatively the control of the population could be done by 
assuming an increase in the mortality rate p, E , but the net result would be 
qualitatively the same. Finally, in the last equation the term 



represents the rate by which infected eggs grow and the term 

Pi (cj — di sin (2-nft + tt)) I&9 (cj — dj sin (2irft + 7r)) , 

as already mentioned, is the fraction of the hatched infected eggs which 
evolves to infectious adults. 

2 An approximated threshold condition 

In the first part of this section we deduce a threshold condition for epidemic. 
The intuition behind the procedures is discussed later on. 

In order to deduce the threshold condition for epidemic we replace the 
non-autonomous system by a autonomous one, by regarding the time on 
the right side of the system ([TJ as a parameter and then carry out a local 
stability analysis. We linearize the second, the fifth, the sixth and eighth 
equations of the autonomous system around a small amount of disease in, 
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Im, im and ^e- 



iLa. — „h±L 



dt **" N 



Pi (cj — dj sin ($)) i E 9 (cj — dj sin ($)) 

If = 9rM(l-^)iM-^ E i E - 

pi (ci — dj sin ($)) (c/ — dj sin ($)) 

where $ = 27r/t + 7r. 

We then examine the stability of the trivial solution of system (J2J), that 
is, i E = 0, Im = 0, %u = and «m = 0, as if the system were autonomous^. 
For this we assume the solutions: 

in = ci exp (At) 
Im = c 2 exp (At) 

(3) 

im = c 3 exp (At) 
i E = c 4 exp (At) 

drop the Heaviside functions by assuming ci > di, and replace (j3J) into 
equation (j2J). The characteristic equation associated to system (J2J) is then 
obtained: 



(2) 
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-(\ + lH + a H + fi H ) a&|g| 



a§f - - ae(~^) x 



Nmft-n) -At ~(^ + / i M) 

ae(-^r) ^gzIil e -Ar o -(A + /x M ) P/ ( Cj -d /S in$) 

^r M (l - ff ; 



-X- He- 
Pi {ci-di sin($)) 

(4) 

If all the roots of equation (jlj) have negative real parts, then the equilib- 
rium without disease is stable, that is, the origin is an atractor.. As shown 
in [3], the first root that crosses the imaginary axis do so through the real 
axis and this happens when 



f>(f\ 2 Nmjt-T!) acxp(-^ M T)k S H (t) 

A W ^ H+aH+flH ) N H (t-n) m m N H (t) 

(5) 



+ 

Vl{ci~di sin*)gr M (i-f^J 
^u(^E+Pli c I- d l sin*)) 



> 1 



Note that the first term in equation (JHJ) is exactly the expression proposed 
in [Hj for the so-called 'basic reproduction number'. 

The intuition behind the above procedure is the following. System (0) has 
'no-mass', that is, it responds to perturbations instantaneously. Therefore, 
we can find the time t at which the stability of the trivial solution of system 
(J2J), that is, %e = 0, in = and %u — becomes unstable. We have numeri- 
cally checked that the time t at which the trivial solution (no-disease) of the 
autonomous system becomes unstable (R > 1) corresponds approximately 
to the moment at which the epidemic takes off, that is, when the epidemic 
in system (P) begins to increase as a result of the introduction of a small 
amount of disease at time t = 0. 
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3 Qualitative analysis of the system's behav- 
iors 



In this section we analyze qualitatively all the possible behaviors of the sys- 
tem when a small amount of disease is introduced into a previously uninfected 
population and when R(t) is in its minimum value (winter time). We do so 
by using R(t) as in equation (jSJ) and by numerically simulating system © 
with parameters values as in table 1. 

Table 1 

We analyze two epidemiological scenarios, one in which R(t), in the ab- 
sence of infection, is most of the time above one, and another in which R(t) 
is most of the time below one. In the first case, if a small amount of infection 
is introduced we observe a pattern shown in figure 1. In the second case a 
small amount of infection introduced generates a pattern shown in figure 2. 

Figure 1 
Figure 2 

In figure 1 the intensity of transmission is relatively low (a = 3.7 days -1 ) 
and we see a first peak followed by a succession of outbreaks forming a 
damped oscillation pattern and the disease disappears. In other words, after 
the first outbreak the infection transmission decreases to levels inferior to 
that of the previous cycle. As the system oscillates and the time interval 
during which R(t) is above the threshold for transmission is insufficient for 
keeping transmission, we have the pattern observed. 

In figure 2 the intensity of transmission is higher (a = 4.3 days -1 ) than 
that shown in figure 1 and, consequently the time interval during which 
the system is above the threshold is larger. In this case, the amplitude of 
the consecutive outbreaks increases until the fraction of immune individuals 
reaches a herd immunity threshold and the disease dies out in a damped 
oscillation pattern. 

We have numerically found that there is an increase in the amplitude of 
consecutive outbreaks with a = 4.3 days -1 , whenever the preceding period 
of time R(t) is above the threshold for transmission is greater than a certain 
time interval (about 190 days). Indeed, with a = 3.7 days -1 , as mentioned 
above, there is a first outbreak and the subsequent peaks formed a damped 
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oscilation and the time interval R(t) is above the threshold for transmission 
is greater than 190 days only for the first peak. 

When we simulate the system with parameters that make R(t) < 1 the 
disease cannot invade the population and disappears exponentially. 

Those are the only three possible qualitative patterns generated by a 
small amount of disease introduced into an entirely susceptible population 
when R(t) is at its minimum value. 

Let us concentrate in the first two patterns, which can better be visualized 
in figures 3 and 4, where the threshold parameter R(t) is shown as function 
of time for each of the above cases. 

Figure 3 
Figure 4 

Other interesting results are shown in figures 5 and 6, in which the time 
oscillation of the 'total amount of disease', d(t), defined as 

d(t) = ^{I H {t)f + {L M (t)f + {I M {t)f + {I E (t)f (6) 

is plotted together with R(t), for both cases of low and high intensities of 
transmission. It can be noted from the figure, that the points in which R(t) 
crosses 1 corresponds, approximately, to maximums and minimums of the 
function d(t). Note also that, in both cases the peaks and troughs of d(t) 
occur slightly after R(t) crosses 1, decreasing and increasing, respectively, as 
if the system has a small 'inertia'. 

Figure 5 
Figure 6 

4 Final comments 

This paper presents a novel, as far as we know, approach to analyze the 
response of a non-autonomous system to a perturbation. This is quantified 
by an approximate expression for the threshold condition that determines 
whether the system will amplify or reduce a small of disease introduced at 
time t. We should warn the reader that we have no mathematical proof of 
the correctness of the threshold expression here deduced but the numerical 
investigation presented above suggests that it is basically correct. 
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A possible application exemplified in this paper is the case of dengue fever, 
in which a seasonal variation in the density of vector mosquitoes determines 
the intensity of transmission. In another paper we used a similar model 
to explain the question of overwintering, that is, how dengue fever survives 
through the winter's dry and cold season. 

Finally, we think that the analysis proposed in this paper could be applied 
to other vector-borne infections and also to some directly transmitted diseases 
that show seasonality in the intensity of transmission. 
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Captions for the Figures 



Figure 1. Number of infected humans for the case of low intensity of 
transmission (a = 3.7 days -1 ). There is a first outbreak resulting from the 
introduction of a small amount of infection in a previously uninfected popu- 
lation followed by a pattern of damped oscillation until the disease disappear. 
Simulation begins at the peak of the winter. 

Figure 2. Number of infected humans for the case of high intensity 
of transmission (a = 4.3 days -1 ). After the first outbreak resulting from 
the introduction of a small amount of infection in a previously uninfected 
population there are subsequent outbreaks with larger amplitudes until herd 
immunity is achieved and the disease gradually disapears. Simulation begins 
also at the peak of the winter. 

Figure 3. The threshold R(t) in the case of low transmission (a = 3.7 
days -1 ). 

Figure 4. The threshold R(t) in the case of high transmission (a = 4.3 
days -1 ). 

Figure 5. The threshold R(t) and 'total amount of disease' d(t) in the 
case of low transmission (a = 3.7 days -1 ). The peaks and troughs of d(t) 
occur slightly after R(t) crosses 1. 

Figure 6. The threshold R(t) and 'total amount of disease' d(t) in the 
case of high transmission (a = 4.3 days -1 ). Again, the peaks and troughs of 
d(t) occur slightly after R(t) crosses 1. 
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Table 1 



Parameter 


Meaning 


Value 


a 


Average Daily biting rate 


see text 


b 


("1 1*1*1*1 i TP i* 

Susceptibility to Infection 


0.1 




TT TV T J 1 T\ jT il*i T"\ i 

Humans Natural Mortality Rate 


4 x 10 _5 days _1 


r H 


TT TV 1 j 1 * t j 

Humans Maltnusian Parameter 


ii 1 

1 days 


k H 


Humans Carrying Capacity 


10 6 


OtH 


T y t t 1 -l r jl*J * TT 

Dengue Induced Mortality m Humans 


1 r\ 1 l 1 

10 Mays 1 


1h 


Humans Recovery Rate 


0.143 days -1 


Ps 


Proportion of non-infected eggs that reach adult phase 


0.15 


cs 


Climatic factor modulating winters and summers 


0.08 


d s 


Climatic factor modulating winters and summers 


0.06 


f 


Frequency of the seasonal cycles 


2.8 x 10~ 3 days" 




atj i j i * j i r *i 

JNatural mortality rate of mosquitoes 


0.263 days 1 


T 


Extrinsic incubation period of dengue 


7 days 


«M 


Dengue induced mortality in mosquitoes 


negligible 




Eggs Malthusian parameter 


50 days -1 


Pi 


Proportion of infected eggs that reach adult phase 


0.15 


Cl 


Climatic factor modulating winters and summers 


0.06 


dj 


Climatic factor modulating winters and summers 


0.06 


9 


Proportion of infected eggs laid by infected females 


see text 


k E 


Eggs Carrying Capacity 


10 6 


He 


Natural mortality rate of eggs 


0.1 days -1 
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